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A spherical  glow  in  low-pressure  nitrogen  for  transition  from  abnormal  to  Townsend  regimes  of  a 
spherical  discharge  has  been  studied.  The  model  includes  non- stationary  drift- diffusion  equations 
for  electrons  and  ions  coupled  with  Poisson  equation  for  radial  electric  field.  The  auto- 
oscillations of  plasma  parameters  are  obtained. 


1 . Introduction 

A spherical  discharge  was  observed  for  the  first  time  in 
1997  [1].  It  was  realized  in  the  gap  between  a central 
point-like  anode  and  surrounding  it  cathode. 

This  paper  presents  the  drift-diffusion  numerical  study 
of  a spherical  glow  discharge  in  nitrogen  in  a wide  range 
of  gas  pressures  and  electric  circuit  parameters. 
Oscillating  regimes  were  obtained  for  transition  regimes 
from  abnormal  to  Townsend  discharge  forms. 


2.Model 

A volume  spherical  discharge  was  studied  in  the  drift- 
diffusion  approximation.  Non- stationary  balance 
equations  for  electrons  and  ions  were  used: 

= |ncUc|Na1()n(E/N)-prccncni , (1) 
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^ +7 d(r2frUi)  = huclNaion(E/N)- Precneni  , 

where  E(r)  is  the  radial  course  of  the  electric  field,  N is 
gas  density,  (3rec  is  two-body  dissociative  recombination 
coefficient  in  e+N2+-^2N  collision,  oc,on  is  ionization 
coefficient,.  Electron  and  ion  fluxes  are  presented  by 
drift  and  diffusion  terms: 

niui  =HiEni-Di^-,  (3) 

The  ionization  and  diffusion  coefficients,  electron  and 
ion  drift  velocities,  Wei=jiejE,  were  taken  in  the  local 
field  approximation.  The  electric  field  was  obtained 
from  the  Poisson  equation: 
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Equations  (1-4)  were  solved  as  the  initial  problem  with 
boundary  conditions  for  electron  and  ion  fluxes  at  the 
cathode  and  anode.  At  the  initial  time  step,  small 
electron  and  ion  densities  (~104cm'3)  were  set  equal  to 
each  other,  and  the  Coulomb’s  distribution  was  set  for 
the  initial  radial  distribution  of  the  electric  field. 

In  the  model,  it  was  assumed  that  a spherical  discharge 
covers  the  entire  cathode  surface  (abnormal  regime  of  a 
glow  discharge,  subnormal  and  Townsend  discharges). 
For  the  numerical  solution  of  Eqs.(l-4),  an  implicit 
finite  difference  scheme  was  used.  As  boundary 
conditions  for  the  Poisson  equation,  there  was  used 
discharge  potential  UdiS,  which  was  recalculated  with  the 
help  of  the  Ohm  law:  Udc=Udjs-RJc,  where  Ups  is  the 


potential  of  power  supply,  R is  the  active  resistance  of 
external  electric  circuit  and  Jc  is  the  conduction  current 
at  the  cathode, 

3.  Results 

Calculations  were  carried  out  for  different  gas  pressures 
(0.1<  p < 3.5  Torr)  and  parameters  of  electric  circuit 
(Ups=  1-3.5  kV,  R=  10-25  kQ).  Radius  of  the  central 
anode  was  fixed  (a=lcm),  the  radius  of  the  spherical 
cathode  was  Rc=12.5  cm.  These  parameters  correspond 
to  conditions  of  [1]  in  which  spherical  stratification  of 
glow  discharge  was  observed. 

As  a result  of  numerical  calculations  the  radial 
distributions  of  electron,  ion,  and  bulk  charge  densities, 
electric  field  and  potential  of  discharge  were  obtained. 


Fig.l.  Relaxation  (1-2)  of  electric  field  at  the  cathode 
(Rc=25cm)  for  different  AUps  = Ups*(p)  - Ups  and  auto- 
oscillation (3)  Decaying  regimes  (curves  4-5)  p=3.5 
Torr,  and  different  power  supply  voltages  Ups. 

It  was  found  that  steady  state  conditions  in  abnormal 
regimes  of  the  spherical  glow  discharges  were  formed  if 
sufficient  voltage  of  power  supply  Ups  was  applied.  In 
this  case,  after  sudden  change  Ups,  decaying  transitional 
oscillation  of  all  plasma  parameters  to  new  stable 
conditions  arose.  In  Fig.l,  the  example  of  such 
relaxation  from  two  different  initial  values  of  power 
supply  is  shown.  The  curves  1 and  2 present  transition 
from  Ups=  1500V  and  Ups= 1350V  to  Ups=1400V, 
correspondently.  It  is  seen  that  the  system  (for 
transitions  from  different  initial  states)  works  out  the 
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steady  state  condition,  which  corresponds  to  the  new 
voltage  of  power  supply. 

For  small  supply  voltage  Ups  (smaller  than  some  critical 
value  Ups ) at  a given  gas  density,  a spherical  discharge 
turns  into  auto-oscillating  (self-maintaining)  regimes  or 
into  decaying  modes.  For  critical  value  Ups , almost 
sinusoidal  oscillation  of  electric  field  and  conduction 
current  are  realized. 

It  should  be  stressed  that  the  critical  voltage  Ups*  is 
realized  only  for  gas  pressures  higher  than  some  critical 
value  p*  depending  on  the  radius  of  spherical  discharge 
Rc.  For  spherical  discharge  gap  Rt,=12.5  cm  and  active 
resistance  of  external  electric  circuit  R=20  kQ,  the 
critical  values  are  p*~1.9  Torr  and  Ups*  = 800  V.  The 
oscillation  frequency  for  this  case  is  f ~ 30  kHz.  For 
higher  gas  pressure  p=3.5  Torr,  Ups  = 1250V,  f = 60 
kHz.  This  regime  is  shown  in  Fig.l  (curve  3). 
Frequencies  of  oscillations  correspond  to  the  time  of 
ions  move  from  the  center  of  cathode  fall  to  the  cathode 
surface  where  they  cause  secondary  emission  of 
electrons. 

In  Fig.2,  radial  distributions  of  electric  field  for 
oscillating  regime  of  a spherical  discharge  at  different 
moments  are  presented  (Ups=1.25kV,  p=3.5  Torr).  It  is 
seen  that  the  amplitude  of  electric  field  changes 
considerably  in  every  point  of  discharge  gap. 


r,  cm 

Fig.2.  Radial  distributions  of  electric  field  at  two 
different  moments  1 and  2 with  interval  At=12,5  ps, 
corresponding  to  anti-phase  of  oscillation. 

It  was  found  that  for  some  regimes  with  Ups  < Ups*(p) 
the  second  oscillating  mode  appeared.  The  frequency  of 
this  mode  is  two  times  smaller  than  the  first  one.  The 
amplitude  of  this  additional  mode  increased  with  the 
increase  of  AUps  = Ups*(p)  - Ups.  The  example  of  such  a 
bifurcation  is  presented  in  Fig.3  for  p=4  Torr,  Ups* 
=1300  V,  AUps  ~ 25  V. 

The  nature  of  such  complex  oscillations  is  connected 
with  the  “breathing”  of  the  width  of  cathode  layer. 

For  larger  AUps  the  oscillating  processes  become 
aperiodic,  and  a discharge  turns  into  the  decaying  mode 
(see,  Fig.l,  curves  4-5). 


Fig.3.  Bifurcation  of  the  period  of  electric  field 
oscillations  Ec(t)  at  the  cathode.  The  curve  is 
superposition  of  two  modes  with  frequencies  ~ 45.5 
kHz,  f2  ~ 22.7  kHz,  and  amplitudes  E(  =52  V/cm,  E2=27 
V/cm. 

4.  Conclusions 

A DC  spherical  glow  discharge  has  been  studied  by 
means  of  self-consistent  drift-diffusion  approximation  in 
a wide  range  of  gas  pressure  and  electric  circuit 
parameters. 

In  transition  regime  between  abnormal  and  Townsend 
regimes  of  discharge  burning,  auto-oscillating  and 
decaying  modes  were  obtained.  All  discharge 
parameters  (conduction  current,  ion  and  electron 
densities,  and  electric  field)  oscillated  in  the  whole 
range  of  the  spherical  discharge  gap  with  the  frequencies 
f~10-100  kHz. 

The  nature  of  such  oscillations  is  connected  with  the 
periodical  ionization  of  molecules  in  the  cathode  sheath 
by  electron  impact,  drift  of  ions  to  the  cathode  surface 
and  secondary  emission  of  electrons  from  the  cathode. 

In  certain  cases,  the  bifurcation  of  the  main  mode  of 
oscillation  was  obtained.  The  second  oscillating  mode 
with  two  time  smaller  frequency  appeared  due  to  the 
doubling  of  the  cathode  layer  width  during  one  half  of 
the  oscillation  period. 

Similar  oscillations  have  been  observed  in  the 
experiments  (see,  for  example  [2])  in  the  transition 
region  from  Townsend  discharge  to  a glow  discharge. 
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